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In this work we present a procedure to infer the mass of progenitors and remnants of Gamma Ray 
Bursts (GRB), starting from the observed energy emitted isotropically and considering the 

associated emission of Gravitational Waves (GW) Eff^ in the different phases. We assume that the 
GW energy of the progenitor Ep^oG is emitted partially during a star collapse, and the residual 
energy is related to the GW energy emitted by the remnant. 

We take a sample of 237 Long GRB, and use an hybrid Montecarlo procedure to explore, for each 
of them, a region of possible solutions of E^^Y s^s a function of the masses, radii, oblateness, rotation 
frequencies of progenitor and remnant and the fraction of energy k emitted as GW by the GRB. 

We discriminate between a Neutron Star (NS) or Black Hole (BH) for the remnant and obtain 
interesting values for the GW emitted by the remnant NS or BH, for the conversion factor k of and 
for the masses and radii of GRB progenitor stars. We also obtain remnant populations with mean 
masses, mean GW frequencies and GRB frequency of GW emission in agreement with the most 
accepted models. 

PACS numbers: 


I. INTRODUCTION 

The leading theories for the origin of long GRB (LGRB) assume that the progenitors are massive stars, such as 
for example Wolf-Rayets (WR) and red giant, collapsing to a compact object (black hole (BH) or neutron star (NS)) 
[b^. Wolf-Rayet stars are candidate progenitors for type Ib/Ic supernovae (SNe), which give birth to either a NS or a 
BH [^. They may be the progenitors of long and soft GRBs, or at least part of them [^. The question is whether 
WRs collapse into a BH, or give birth to a type Ib or type Ic SN. This is relevant, because in this paper we aim to 
classify the GRB observations according to their progenitor mass. In the scenario of possible GRB progenitors, type 
Ib, Ic and Ib/c supernovae (SNe) or hypernovae are produced [IMllIl- It is realistic that core-collapse SNe induce 
LGRB signals when stellar cores undergo fast rotations before the collapse. The bulk of matter then forms a torus 
around the newly-formed central compact object and, following the accretion of the torus, feeds the GRB [13, 

As it is well known, the first non-electromagnetic signal that can be related by a GRB is a gravitational radiation 
emission in the form of GWs. Such perturbations are not generated by the burst itself, but by the progenitor through 
the central engine [s^- Furthermore, in the case of LGRBs, the GW signal should arrive soon before the burst, i.e. 
during the core-collapse phase or during the production of the central BH which could be originated also by merging 
processes. For the SNe-GRB connection, several papers have been written; a comprehensive discussion on the problem 
can be found in (b^ and references therein. Considering progenitors, particularly in the case of LGRBs, there are 
no final models for the GWs emission but several remarkable papers have been published pointing out the problem 
@1 [H, m, in, 113,113, [bll ■ The expected signal would lie in the range 10^ -L 10^ Hz and, in principle, could be detected 
only for events up to a few tens of Mpc (13|. The limitation in distance is the main and critical problem for the 
detection of possible GW signals related to GRB progenitors. In this paper we suppose that the progenitor masses, 
at the origin of LGRBs, could be related to the compact objects resulting as remnants of SNe or collapsar processes. 
The resulting remnant can be either a NS or BH objects taking into account the emitted energies as GW in the GRB 
prompt and in the afterglow phase. The GW production can give hints on the total mass of the collapsing system 
by estimating the global emitted gravitational energy. By relating the energy emitted in GW by the progenitor to a 
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fraction k of the electromagnetic GRB prompt emission energy, and the energy emitted in GW by the remnant, one 
could trace the progenitor mass. 

The main steps that we will follow are: 

• We assume that the energy released by the GRB progenitor during the collapse is partially converted into GWs, 
i.e. we assume that the order of magnitude of gravitational energy is a certain fraction k of the energy emitted 
in the electromagnetic 

• As a second step, we constrain the progenitor masses assuming that the energy balance of progenitor in the 
collapsar phase will be: Ep^QQ = E^^^ + Ep^pf, where Ep^QQ is the energy emitted by the progenitor as 
GW’s, Efg^ is the energy converted in GW during the GRB prompt, and finally Ep^pp is the GW energy 
emitted by the remnant. We are aware that the energy emitted in GW by the remnant is not in principle 
directly related to the energy of the progenitor, but this is a rough estimate to correct the total GW emission 
of the GRB. On the other hand it can be roughly justified by the consideration that some of the energy of the 
progenitor goes into the kinetic energy of the remnant, e.g. in rotational energy, and thus in the latter GW 
emission. 

• We compare the dimensionless amplitude ho of GWs, during the collapsar phase, with the current sensitivity 
curves of the ground interferometric GW antennas, like LIGO and VIRGO, and the future interferometric GW 
antennas, like Advanced LIGO, Advanced VIRGO and the forthcoming Einstein underground interferometric 
GW antenna (ET) (see d, H, [^); 

In this context we choose progenitor masses in the interval IOMq < M < 60Mq, considering the more massive 
progenitors that can also be 60Mq. Therefore, in this scenario, we take also into account the Wolf-Rayets with 
C/O < 1 (WO) stars, that are produced in a limited mass range (around 60Mq) and only at low metallicity 
{Z = 0.010). In the results presented in they end up with a NS or a BH as a remnant depending on the initial 
progenitor star mass, which is between 15Mq < M < 32Mq in the first case, while in the latter 32Mq < M < 
60Mq. However, we here point out that, according to [s^, a scenario of 35 solar masses is preferred. Moreover, 
among the progenitor types they select in the mass ranges between 12M0 < M < 20Mq Red Super Giant (RSG), 
while for the other cases WC and WNL. present grids of stellar evolution models at Z = 0.004,0.002,0.001, 
and 0.00001, for rotating magnetized stars in the initial mass range 12 Mq < Minu < 60Mq, and with initial 
equatorial rotation velocities between zero and 80% of the Keplerian value (0 < < 0.8). For Z = 0.001, 

the mass range of GRB production is 12 Mq < Minn < SOMq, and 20 Mq < Mmn < 4OM0. In addition, 
GRB progenitors have a more massive helium envelope for slower semi-convection, on average. At Z = O.OOOOl, 
on the other hand, the rapidly rotating stars with M > QOMq form CO cores of M > 4OM0, which may be 
unstable due to the pair instability (cf. js^l)- The precise CO core mass limit for the pair instability may 
increase with higher core angular momentum [^. We take these results as a reference for our study to achieve 
probable dimensionless GW amplitude and the GW emission frequency connected to the GRB. In fact, the mass 
M fv 25M0 is near the transition from NS to BH formation [s^, the exact boundary depending on rotation and 
mass loss [^, while we, conservatively, chose the star of M Ri IOM 0 [^ . 

The paper is organized as follows. In Sec. 2, we discuss the possible GW emission associated to LGRBs. The goal 
is to show that the energy emitted in GW by the collapsing star can lead to the determination of the progenitor mass 
and to a suitable estimate of the fraction /c << 1 of the electromagnetic energy E^^^ that is emitted in GWs, and 
the energy emitted in GWs by the remnant. In Sec. 3 we discuss the algorithm we used, to obtain the progenitor and 
remnant parameters from the observed Ep,^^ ± . We used the Ep.^^ values observed by Swift , which span 

over a very wide redshifts range. Results of calculations are reported in Sec. 4 while Sec. 5 is devoted to discussion 
and conclusions. In particular we discuss the initial conditions required to obtain collapsars and LGRB. 


II. GRAVITATIONAL WAVES EMISSION ASSOCIATED WITH LONG GAMMA RAY BURSTS 

The progenitors of LGRBs can be considered as the final step of a stellar collapse and their modeling is rather 
complicated since cumbersome three-dimensional simulations are required. Many simulations have been performed 
for intermediate cases of core-collapse of SNae, which should provide GW bursts [^. Gravitational collapse leading to 
the formation of a NS has long been considered an observable sonrce of GWs. During the core collapse, an initial signal 
is expected to be emitted due to the changing of mass distribution. A second part of the GW signal is produced when 
gravitational collapse is halted by the stiffening of the equation of state above nuclear densities and the core bounces, 
driving an ontward moving shock, with the rapidly rotating proto-neutron star (PNS) oscillating in its axisymmetric 
normal modes. In this paper we show that, at least for our GRB sample, taking into acconnt the GRB distances 
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(« Gpc, well beyond the horizon of GW antennas), with the Advanced VIRGO and LIGO we will not be able to 
see GRBs within a certain ranges of frequencies and amplitudes, also according to the energy conversion factor from 
electromagnetic to GW that we estimate. 

As a first step, we start by evaluating the dimensionless GWs amplitude h^. We remember that, in cases where the 
GWs emission comes from a core-collapse process, we cannot predict the exact shape of the emitted GWs signal, but 
we can have an estimate starting from the energy flux. Hereafter we will use the following notation: 

poo 

h{t)= (1) 

J oo 

The energy flux of GWs, defined as the power emitted per unit of surface at a large distance from the source is [^ : 


•^GW = 


d 2 TT’ZSO 

^GW 

dtdS 


IGttG 




( 2 ) 


where <...> indicates a temporal average over a large enough number of periods. The total energy emitted assuming 
isotropic emission is then 


= An J dt Egw, 


( 3 ) 


where D is the distance to the source. We can write the power emitted in GW by the GRB (the gravitational 
luminosity), by inserting eq.([T]) in eq. ([3|) to obtain: 


pGW 

^iso 


n^c^ 

AG 




fho{f)e 


— 27Tift 


df 




( 4 ) 


being 


hUf) = + (5) 

If the frequency distribution of the GRB is peaked around a frequency few, as e.g. in a sine-Gaussian waveform, eq. 
@ can be written as: 


= ^D^2fawho{tr. ( 6 ) 

Considering that the only observable in the GRB emission is the flux of electromagnetic energy Efl^^ (in a certain 
band), we make the assumption that the energy emitted in gravitational waves is a fraction k of . So we have: 


pGW 

^iso 


kEgr 

Tqo 


where Tgg is the time in which the 90% of E^^^ is emitted, and so: 


( 7 ) 


ppGRB 

The fraction k will also be a tunable parameter useful in order to compare the present and foreseen GW inter¬ 
ferometric antenna sensitivities {i.e. Adv. VIRGO, Adv. LIGO and the forthcoming underground GW Einstein 
Telescope) to the mean scaled GW amplitudes. 

Let us consider a progenitor star of mass Mppc, radius Rprg and oblateness eppoGi rotating with a frequency 
fpRG = O-^ fpRGn-u,, and a NS as a remnant after the process of GRB emission with parameters M^s, Rns^ ^ns and 
Jns- We have that [3^ : 


h 


PRG 

GW 


An'^G ^s^^fpRGcw 


D 


iPROG 


( 9 ) 


h 


NS 

GW 


An^G 


' tNS f2 
^3 JgWns 

D 



( 10 ) 
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The momentum of inertia for a spherical rigid object (considering negligible the quadrupolar oblateness e) is given 


by: 

^3^ = -^MnsRns^ 

(11) 


= -MppgRprG ■ 

(12) 


The hypothesis of a rigid body also for the progenitors, can be accepted if one considers that magnetic fields can 
enhance their stiffness. Taking into account the total energy balance of the process, we can roughly say that the 
energy emitted in GW by the GRB is given by the difference of the energies emitted in GRB by the progenitor and 
the remnant resulting in a similar relationship between the squares of GW amplitudes. From eqs. (ED, (ED, (UnD, (ED 
and (|12L we can very roughly write, in terms of the only observables, i.e. and T90: 


jpGRB 

^iso 




25c^k 


so we get a function useful to solve our problem, that is 

^Gw = fik,ei\rs,Mi\rs,RNS, fNS,epRG,MppG,RpRG, fpRG, /grb) ■ 


(13) 


(14) 


To evaluate the physical parameters involved in the GRB process to the remnant phase, we should estimate the 
= 10 parameters in eqs. (fT^ having only one observed data i.e. . However the assumption we make for the 

progenitor models are very well posed in literature so that we can set very reasonable ranges for 9 parameters out of 
10, therefore the only key parameter we infer is the scale factor k. 

In the next section we will put the bases to solve this problem using a key idea and a suitable algorithm. 


III. GENETIC CONTROLLED RANDOM SEARCH ALGORITHM FOR GRB PROGENITOR MASSES 

ESTIMATE 

Our aim, as stated above, is to find the best progenitor mass starting from the observed E^^^. Given eq. (|13D we 
see that, from the observed Ef^^^ we must infer ten parameters: 

1. The factor k to evaluate the scaled E^^ from the observed E^^^; 

2. The emission frequency Jgw^^ 

3. The estimate of the progenitor quadrupolar oblateness. According to interferometric measures and stellar 
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evolution grids, it can be put in the interval 0.01 - being the upper limit obtained from the star Roche lobe 

limit, at the critical velocity (i^ . 

4. The mass of the progenitor. According to Grids of stellar models with and without rotation (^ . we assume 
the progenitor mass in the interval 10 ^ 60Mq. 

5. The frequency of GW emission of the progenitor. According to the Grids of stellar models (^. we assumed 
it in the range: 10“^° ^ 10~^Hz inferred from the critical velocity according to the Roche model (^ . 

6. The estimates of the remnant quadrupolar oblateness. According to the current NS model and the corresponding 

Equation of State (EoS) it can be set in the interval 10“^ 10“® (i^ . 

7. The mass of the remnant. According to different EoS models the mass of the remnant could be in the range 
1.0 ^3M0 m,!!! or for a stellar BH in the range 1.0 IOMq [gOI ■ 

8. The radius of the progenitor. According to the Grids of stellar models and references therein) we can 

assume stellar radii of the progenitor in the interval 0.9 ^ 140i?Q from stellar model of mass ranges 10 ^ 60Mq. 

9. The radius of the remnant. According to different EoS models, the radius of the remnant could be in the range 
bKm < R < 15Km, [i^ or for a stellar BH in the range 5Km < R < 30Km. 

10. The frequency of GW emission of the remnant. From the known frequencies of electromagnetic energy emission 
by pulsars we can consider GW frequencies Jns in the range O.IO ^ lOOOH^ 
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Table I: Model parameters ranges used according stellar evolution model of 


Name 

MpROG [Mq) 

Rprog {Rq) 

Muemn (Mq) 

RRemn 

f Rerun (^■^) 

Model 1 

10t60 

0.9 T 140 

1t3 

6t 15 

0.1 T 1000 

Model 2 

10t60 

0.9 T 140 

It 10 

6t30 

0.1 T 1000 

Model 3 

IOtSS 

0.9-;-25 

1t3 

6t 15 

0.1 T 1000 

Model 4 

10t35 

0.9-;-25 

It 10 

6t30 

0.1 T 1000 


Table II: Parameters common to the different model of Tab. |I] useful to initialize the domain of possible solutions. 


Model 

k 

fclpB (Hz) 

bJpRoc {rad/s) 

epROG 

^Remn 

Model C 

10~® T 10"'‘ 

50 T 800 

10"'^ T 10“^“ 

10"^ T 1 

10“® T lO""' 


To solve the problem we use an algorithm that is capable to find the global minimum of a multivariable function. We 
used a simplihed simplex algorithm: the controlled random search (CRS) algorithm and, to improve its performances, 
we used a genetic modification of the search procedure making the software more resilient to local minima 
). In Fig. [T]we show the flux diagram of the algorithm. 

First of all we must obtain an objective function to minimize, using the reduced Xobs Nqrb — N degrees of 
freedom, where = 10 is the number of free parameters: 


^Ob Xobs 


1 


Ngrb 

Ngrb — N ^ 
1—1 


{E,, 


— E,,, 


(15) 


where is given by Eq. (11311 . It worths noticing that we have only one observed datum i.e. Eisooba EuEiaa 

GRB, in other words one equation and ten unknown parameters to estimate. We solve the problem by a Montecarlo 
procedure, for each observed GRB, we generate a large number Ngrb of Eiso values extracted randomly from a 
Gaussian distribution with mean value Eisooba standard deviation crEi^a given by the observed measurement 
error: 


Eisogar. = E^sooba + Ei8orandn{NG rb) , (16) 

then we evaluate the reduced x%t — Xobs with Ngrb — N degrees of freedom at 5% of conhdence level {xju ^ 0.86). 
Finally we evaluate the mean Eisog^nM Eisoga„ lo check the relative difference between the Observed (O) and 

calculated (G) values: O — C = of the four regions of possible solutions, hereafter indicated as model 

1...4, with parameters in the interval shown in TablD The parameters common to the different models of TablH to 
initialize the domain are shown in Tab. im 

Using these ranges of parameters we analyzed a sample of 237 GRB and the results are treated in the next section. 
In Fig. [T]we show the flux diagram of the algorithm, and an example of convergence curve using the value of the Xobs 
as a function of the number of iterations. The red circle, represents xju ^ Xobs ~ giving a confidence level of 
5% for the test case of GRB 060218 that we will analyze later. 


IV. GAMMA RAY BURSTS SAMPLE AND RESULTS ON GRAVITATIONAL WAVE EMISSION 

We selected a sample of 237 long GRBs with known redshifts that we then analyzed, detected by the Swift satellite 
[23 from January 2005 to May 2014. This sample is an extended sample presented in [l4 - [l^ where GRBs with plateaus 
are presented. We use the GRB Coordinates Network 0 for redshifts and exclude GRBs with non-spectroscopic 
redshifts. In Tab. m the sample GRB data range, that we analyzed, is shown. 

Considering the above mentioned sample of long GRBs and according to the hypotheses we made, it is possible to 
infer the physical parameters of the GRB progenitor with their errors. The mean results of our analysis are shown in 


Table III: Sample GRB data range used in this work 


{erg) 

a^GRB {erg) 

T90(s) 

orgo (s) 

redshift 2 

D X 10^® (cm) 

2.5 X 10^® T 1.1 X 10®® 

3.9 X 10'‘’' T 1.7 X 10®'‘ 

2.2 T 844 

0.003 T 3 

0.014 T 8.2 

1.8 T 2500 
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Figure 1: On the top panel we show the flux diagram of the algorithm, and an example of convergence at confidence level of 
5% Xfit ^ Xobs ~ 0.86 (red circle) for the test case of GRB060218. On the bottom the four region of possible solutions for the 
progenitor masses and radii are shown according to the used parameters range. 


Tab. IIVI and Tab. |V] whilst in Fig. [5] the residuals O — C% of the fitting against the EisoQ^^ are shown in log scale for 
the four models with parameters interval shown in TabUand Tab. In] It is evident that the agreement is fairly good. 
The mean results of parameters common to the different models of Tab llVI are shown in Tab. El 
So we have the astrophysical parameters of the possible progenitors and remnants and the GW emission frequency 
of GRB. The scheme of our procedure is shown in Fig. [3] 

The results obtained so far can be analysed in order: 

• to discriminate between stability or instability condition for the progenitor, for the region of possible solutions 
spanned by genetic-price algorithm; 

• to discriminate between NS or BH solution for the remnant; 

• to discriminate among the possible region of solutions obtained according to the models in analysis. 
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Table IV: Mean results for model parameters obtained from genetic-Price algorithm 


Model 

MpROG (Mq) 

Rprog [Rq] 

MRerun {Mq) 

RRemn 

fRemn 

Model 1 

41 ±5 

50 ±12 

2.3 ±0.7 

10 ±2 

504 ± 260 

Model 2 

34 ±3 

36 ±11 

5±2 

16 ±6 

782 ± 300 

Model 3 

25 ±4 

11 ±4 

2.0 ±0.5 

9±3 

495 ± 200 

Model 4 

25 ±4 

10 ±4 

5±2 

15 ±6 

493 ± 230 


Table V: Mean result Parameters common to the different models of Tab. |IV 


Model 

k 

fa'RB {Hz) 

OJPROG {rad/s) 

^PROG 

^Remn 

Model 1 

2± 1 X 10"® 

602 ± 255 

3± 2 X 10"® 

0.15 ±0.06 

3± 2 X 10"® 

Model 2 

3±2 X 10"® 

593 ± 264 

5±2 X 10"® 

0.21 ±0.09 

2± 1 X 10"® 

Model 3 

4± 2 X 10"® 

628 ± 270 

8± 5 X 10"® 

0.23 ± 0.09 

1±0.8 X 10"® 

Model 4 

4± 2 X 10"® 

618 ± 278 

9±6 X 10"® 

0.22 ± 0.08 

5±3 X 10"® 


a. - Stability or instability condition for the progenitor. To discriminate between stability or insta¬ 
bility condition for the progenitor, within the region of possible solutions spanned by genetic-price algorithm we 
used the ratio between the progenitor star angular velocity and the Keplerian velocity (assumed as critical veloc¬ 


ity) UIpRG < ^crit = 


GM 


We compared the unity ratio i.e. Pcrit = 1 and the escape ratio pesc = 


Pcrit 


to the 


ratio between the actual angular velocity of the progenitor and the corresponding critical one. In Figl31 the ratios 
Pcrit = —are shown to ascertain whether the progenitor star is rotating at an angular velocity greater than the 


Keplerian velocity. The results are shown, for the model 3 on the left and the global ones on the right with the error 

^PRG 

bar on pcrit- The black line represents pcrit = 1 whilst the magenta line is the escape ratio Pescape = — 7 =-■ As it 

Pcrit 

is possible to see we have values below the critical and near the escape ratio. It is evident that the dynamical state 
of the progenitors is essentially the same for all the models. The red and the green squares are the p of GRB030329 
and GRB060218 respectively, taken as test cases. 

b. - How to discriminate between NS or BH solution for the remnant. It is worth noticing that the 
remnant, accor ding to the different models, could be either a millisecond pulsar or a rotating BH as it is foreseen from 
models f|4H.lOT.I^ and references therein). To ascertain whether the remnant is a BH or a newborn millisecond pulsar 

h 

like NS (see Fig. O we used the gravitational parameter psch =--— the ratio between the Schwarzschild 

Rr emn 

^GM 

radius, Rsch = — 5 —> and the remnant radius. We used the Schwarzschild radius for the comparison because the 

c . 

difference with the inner Kerr ergosphere radius is surely within the errors in the determination of remnant mass, 
angular velocity and radius of the remnant. In Fig. [5] the black and the blue circles are the results of Model 2 and 
Model 4 (BH zone) whilst the magenta and red circles are the results of Model 1 and Model 3 (NS zone). The different 
symbols (see legenda) are the psch of GRB030329 and GRB060218 respectively. The green line is the limit between 
the BH {psch > 1) and the NS {psch < 1) zone. It is evident that we have two sets of different possible remnants NS 
or BH. 

We are also aware that, in our procedure, we treat the emission of GW from a BH in the same way as the emission 
from a rigid body. Obviously this is not the case, but as far as we are interested in the transition from one state to 
the other, the assumption can be considered as a limiting case. 

c. - How to discriminate among the possible regions of solutions according to the models. In order 
to probe the detectability of the hypothesized GW signals from progenitors, GRB and remnants, we compare the GW 
amplitudes obtained for the different models with the sensitivity curves of present and future GW interferometric 


Table VI: Sample GRB data for test on GRB030329 and GRB060218 


Name 

{erg) 

o-^GRB {erg) 

T90s 

ergo (s) 

redshift ^ 

lum. dist. D {cm) 

GRB030329 

1.80 X 10®^ 

7.00 X 10®° 

41.02 

0.02 

0.1685 

1.81 X 10^^ 

GRB060218 

2.54 X 10"^® 

5.22 X 10'*'^ 

128 

0.2 

0.0331 

4.42 X 10°® 










































Figure 2: Residuals O — C = fitting against Log{Eisoobs) 

Eisoobs 

shown in Tabl!] It is evident that the agreement is rather good. 


for the four model with parameters interval 


antennas: VIRGO, LIGO and Advanced LIGO, Advanced VIRGO and Einstein Telescope. We also used those 
comparisons to try to discriminate among the different models according to the detectability of the GRB and remnants 
of our sample. 

In Fig. [6] on the top and bottom panels we show the GWs amplitudes we obtained for Model I and Model 3 from 
GRB sample (blue and magenta circles respectively) overlapped to VIRGO and LIGO sensitivity curves assuming the 
computed kc, i- e. energy conversion of into GW energy. The GW amplitudes of GRB0602I8 and GRB030329 

is also reported on the plots. The inspection of the figure shows that the Eq^ computed amplitudes, could be observed 
only with the future facilities specially the Einstein Telescope, some being at the limit of the foreseen sensitivity for 
Advanced LIGO, Advanced VIRGO. Specifically, VIRGO and LIGO performances are shown respectively from the 
scientific run (VIRGO) VSR2 and the scientific S'G run (LIGO), compared with the respective target curves [U, Hj. It 
is worth noticing that the amplitudes of GRBqw whatever the Model are practically in the same frequency ranges. In 
the figures there is an interesting prediction of a remnant millisecond population of either of millisecond pulsar (Model 
1 Model 3: red and green circles top) or BH (Model 2 and Model 4: red and green circles bottom). So, according 
to our approach, we inferred the progenitor mass according to the hypothesis of generation of either a NS or a BH 
after GRB explosion from a progenitor star, whose type we may infer from evolution star models, using the physical 
parameters we computed. From Fig. [Slit is easy to see that in some cases we can discard the results because, if the 
model were true we should have already observed GW from the remnants with the present sensitivities of VIRGO 
and LIGO so we can discard Model 1 Model 2 and Model 4. We can affirm so that, according to the remaining Model 
3, the remnant must be NS and the progenitor masses and radii are in agreement with the hypothesis that the GRB 
progenitor must be WR stars. Last but not the least LIGO and VIRGO are presently undergoing to an upgrade to 
significantly increase their sensitivity, so, the chances to detect a coincident emission of GWs with a GRBs will be 
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Figure 3: Astrophysical parameters determined with our hybrid Montecarlo procedure and scheme followed to analyze the 
results. 



Figure 4: The ratios pcrit = shown to ascertain whether the progenitor star is rotating at an angular velocity greater 

than the keplerian velocity. The results are shown, for the model 3 on the left and the global one on the right with the error 
bar on pcrit- The red and green squares are the p s of GRB030329 and GRB060218 respectively. The black line marks pcrit = 1 
whilst the magenta line is the escape ratio pescape = As it is possible to see the dynamical state of the progenitors is 

V 2pcrit 

practically the same for all the models. 
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Figure 5: We show here the gravitational parameter psch = 


Rsch 


being Rsch = 


2GM 


the Schwarzschild radius and 


RRemn 

the remnant radius in order to discriminate whether the remnant is a BH or a newborn millisecond pulsar. The black and the 
blue circles are the results of Model 2 and Model 4 (BH zone) whilst the magenta and red circles are the results of Model 1 and 
Model 3 (NS zone). The different symbols (see legenda) are the psch of GRB030329 and GRB060218 respectively. The green 
line is the limit between the black hole {psch > 1) and the NS {psch < 1) zone. It is evident that we got two sets of different 
possible remnants i.e. NS or BH. 


substantially improved and the prevision of a millisecond remnant population (NS), could be experimentally tested 
also in connection with the estimated conversion energy parameter kc- 


V. DISCUSSION AND CONCLUSIONS 

We have considered the hypothesis that the energy emitted as gravitational radiation from a core-collapsing system 
can be obtained by scaling EQ'^g by a scale factor k ranging from 10“® to 10“® thus admitting that there has been a 
conversion of mass into gravitational energy. From this assumption, it is possible to estimate the progenitor masses. 

From our GRB sample, we analyze the case of GRB060218 whose progenitor could be a WR star Q and the 
case of GRB030329 and references therein. We consider these GRB as test cases, because they could constitute 
an independent support for the newborn magiietars scenario, proposed to account for shallow decays or plateaus 
[T3. . The above mentioned support comes from the observation of SN2006aj, associated with the nearby 

sub-energetic GRB060218, suggesting that the supernova-GRB connection may extend to a much broader range of 
stellar masses than previously thought, possibly involving two different mechanisms: a ’collapsar’ for the more massive 
stars ( 39 I - I 41 I I and a magnetars scenario. In Table IVTIl for GRB060218 and GRB030329, we show, with their errors, the 
progenitor masses in Mq, progenitor radii Rprog bi Rq, remnant masses, radii and GW remnant emission frequency 
of both GRB030329 and GRB060218 respectively, as inferred from our procedure. In Table IVlin we show the k of 
energy conversion, progenitor oblateness epRCi progenitor GW frequency lurrg {rad/ s), remnant oblateness CRemm 
of both GRB030329 and GRB060218 according to the different models. For the case of GRB060218 it is easy to see 
that both the progenitor masses and radii computed according to model 3 and 4, are in good agreement with the 
estimations by Q who states that observations provide strong evidence that the GRB progenitor was a WR star, 
being the star radius definitely smaller than 5 x lO^^cm, i.e. smaller than the radius of the progenitors of type II 
SNe, like blue supergiants (4 x 10^^ cm for SN1987A or red supergiants (3 x 10^^ cm). In Fig. [SlGW amplitudes hiso 
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Figure 6: On top panel GW amplitudes for the GRB sample are shown for models 1 and 3 (blue and magenta circles 
respectively) together with GW emission by remnants (NS for these models, red and green circles) and overlapped to VIRGO 
and LIGO sensitivity curves assuming the computed kc, i. e. energy conversion of into GW energy. The GW amplitudes 

of GRB060218 and GRB030329 are also shown. On bottom panel the same is shown for models 2 and 4 (BH remnants) VIRGO 
and LIGO performances are shown respectively from the scientihc run (VIRGO) VSR2 and the scientific 5'6 run (LIGO), 
compared with the respective target curves P, Ql- An interesting prediction of a remnant millisecond population of either of 
millisecond pulsar {Model 1 — Model 3: green and red circles top) or BH {Model 2 — Model 4: red and green circles bottom) 
can be noticed. 


(see legend) of the remnants of GRB060218 according to the different models are placed over the sensitivity curves 
of GW interferometric antennas in order to discriminate among the different models on the base of the detectability 
of the GRB and remnants. This allows to discard all models except Model 3. Finally we want underline that we are 
proposing a new method of analysis of GRB progenitor mass and, if the redshift is known from Eiso it is possible to 
infer quickly the mass either of the progenitor or the remnant. 
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Table VII: Progenitor mass in Mq, progenitor radius RproGq- Remnant masses, radii and GW remnant emission frequency 
of GRB030329 and GRB060218 according to the different models 


GRB 030329 

zbcr 

MpROG 
±10 Mq 

Rprog 

±5Rq 

MRemn 

±0.8 M© 

RRemn 

±4 (km) 

fRemn 

±2S0 Hz 

Modello 1 

39 

23 

2 

10 

507 

Modello 2 

38 

19 

5 

11 

495 

Modello 3 

24 

10 

2 

10 

491 

Modello 4 

24 

9 

7 

12 

496 

GRB 060218 

±17 

MpROG 
i8 Mq 

Rprog 
zb5 Rq 

MrmN 
±0.5 Mq 

Rrmn 
±3 km 

/rmn 
±250 Hz 

Model 1 

27 

49 

2 

10 

514 

Model 2 

24 

18 

5 

15 

459 

Model 3 

21 

8 

2 

8 

393 

Model 4 

21 

8 

5 

15 

460 


Table VIII: k of energy conversion, progenitor oblateness epRc, progenitor GW frequency luprc {rad/s), remnant oblateness 
^Remn , of GRB030329 and GRB060218 according to the different models 


GRB030329 

k 

pGW 

Jgrb 

OJPROG 

epRG 

^Remn 

±cr 

±2.6 X 10“® 

±270 Hz 

±4 X 10“® rad/s 

±0.07 

±3 X 10"'^ 

Model 1 

5 X 10“® 

542 

5 X 10"® 

0.14 

5 X 10"'^ 

Model 2 

5 X 10“® 

536 

5 X 10"® 

0.13 

5 X 10"'^ 

Model 3 

5 X 10“® 

596 

8 X 10"® 

0.15 

5 X 10"'^ 

Model 4 

5 X 10“® 

613 

8 X 10"® 

0.15 

5 X 10"'^ 

GRB060218 

k 

pGW 

Jgrb 

OJPROG 

epRG 

^Remn 

±17 

±2.5 X 10“® 

±260 Hz 

±1 X 10"® rad/s 

±0.08 

±3 X 10"'^ 

Modello 1 

5 X 10“® 

387 

1.4 X 10"® 

0.10 

7 X 10"'^ 

Modello 2 

4 X 10“® 

400 

2.4 X 10"® 

0.13 

2 X 10"® 

Modello 3 

56 X 10"® 

383 

3.1 X 10"® 

0.12 

2 X 10"® 

Modello 4 

6 X 10“® 

427 

2.7 X 10"® 

0.12 

8 X 10"'^ 
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